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We interpret recent measurements of the zero field muon relaxation rate in the frustrated mag- 
netic pyrochlore Dy2Ti20y as resulting from the quantum diffusion of muons in the substance. In 
this scenario, the plateau observed at low temperature (< 7 K) in the relaxation rate is due to 
coherent tunneling of the muons through a spatially disordered spin state and not to any mag- 
netic fluctuations persisting at low temperature. Two further regimes either side of a maximum 

relaxation rate at T* = 50 K correspond to a crossover between tunnelling and incoherent acti- 
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vated hopping motion of the muon. Our fit of the experimental data is compared with the case of 
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The recent measurement of the zero field /iSR relaxation in Dj^T^Oy spin ice (DTO) 
by Dunsiger et al. [Ij joins a long series of puzzling experiments on a diverse range of 
frustrated magnetic materials over roughly the last fifteen years [IHH]. They indicate a 
relaxation of the spin asymmetry of the muons after they are implanted in the sample in 
a spin-polarised state. This persists down to the lowest observed temperatures with little 
temperature dependence in the relaxation rate below a temperature varying from 1 — 10 K 
depending on the compound. 

The interpretation of this relaxation has long been a matter of discussion. At its center 
lies the question whether the origin of the dynamics in each case is due to intrinsic magnetic 
fluctuations in the material, or whether the implanted muons are instead more than merely 
passive probes of the magnetism. 

The latter is a realistic possibility as the muon couples not only weakly to the magnetic 
degrees of freedom via its spin but also potentially much more strongly to electric degrees of 
freedom via its positive charge. As we argue here, this can give rise to new physics interesting 
on its own right, which is in turn elegantly probed via the magnetic degree of freedom. 

In DTO, previous //SR measurements were made pH] and an important debate [TOfTS"] 
concerning the origin of the muon spin relaxation in spin ice is developing. In spin ices, the 
moments have an Ising anisotropy and the interactions are frustrated leading to the onset of a 
highly degenerate spin ice state at low temperatures that is signalled by a heat capacity peak 
at around 1 K. The dynamics in this material has been explored using several probes besides 
/iSR [HE] including susceptibility [ToTHS] neutron scattering [T9"H2"T] . magnetocaloric effect 
[22], magnetization relaxation experiments [18] . nuclear forward scattering [23], NQR [21] 
and three distinct dynamical regimes have been found. Above around 15 K, the dynamics 
follows an Arrhenius law controlled by a gap to excited crystal field levels of several hundred 
Kelvin [20]. Between about 1 K and 15 K, the dynamics is dominated by tunnelling between 
magnetic configurations and the temperature dependence is correspondingly weaker than at 
higher temperatures [25]. Below about 1 K, the timescales greatly increase, but there is some 
evidence for a second Arrhenius regime in a.c. susceptibility [17] although the moments are 
static on neutron timescales |21j . Between 2-16 K, large timescales (^ 10~ 3 s) have been 
observed by a.c susceptibility measurements [15], and between 30-90 K, nuclear forward 
scattering of synchrotron radiation experiments [23] give characteristic fluctuation times 
between 10~ 7 and 10~ 10 s between 30 and 70 K respectively. 
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FIG. 1: Sketch of the temperature variation of the hopping rate of the polaronic muon within the 
quantum diffusion scenario. 



The existence of several dynamical regimes appears at odds with the almost featureless 
/xSR relaxation 1/Ti below 7 K observed in Ref. pQ in a zero-field experiment. Moreover, 
since the moments are, to an excellent approximation, Ising-like there should be insufficient 
spectral weight at low temperatures to bring about muon spin relaxation. 

In our scenario we assume that the magnetic spins are to a first approximation frozen 
at low temperatures T < 70 K relative to characteristic times experienced by the muons 
and that the relaxation mechanism of the muon is due to its diffusion through the static 
disordered magnetic background |5Uj. So far, such an hypothesis has not been considered 
in insulating oxides. On the contrary, muons are generally believed to be localized in such 
compounds. However, the diffusion of muons has been observed in many substances includ- 
ing some metals and muonium (a bound state of an electron and a muon, denoted by Mu) 
diffusion has been seen in some ionic insulators including KC1, NaCl and GaAs ( [2TH30] ; see 
[31] for a review). The fact that Mu seems not to have been observed in magnetic insulating 
oxides such as Dy2Ti2C>7 is an indication that its formation is screened by the dielectric 
constant of the insulator (e s ~ 65 and e^ ~ 5 [32], where e s and eoo are respectively the 
static and high frequency dielectric constants). Such a high value of the static dielectric 
constant also indicates that the optical phonons have important effects on the muon: the 
muon-phonon interactions with both optical and acoustic phonons almost certainly lead to 
the formation of muonic polarons [21] ■ 

In the following, we assume the existence of muonic polarons in Dy2Ti2C>7 and consider 
the possibility that these polarons delocalize at low temperatures. In particular, we believe 
that the diffusion of just such a polaron offers a scenario (Fig.l) within which the data of 



Dunsiger et al. [1] may be interpreted. 

The remainder of this paper is organized as follows. First, we expose how the muon spin 
relfects the muon diffusion through the disordered spin ice state. Next, we discuss the differ- 
ent regimes of polaron motion, together with a fit to the zero field /iSR data on Dy^T^CV 
The fitting parameters are then compared to independently obtained estimates and found 
to be in satisfactory agreement. We close with observations of how, and to what extent, 
muon diffusion should show up in related compounds. 

As previously mentioned, motion of polarons, based on muonium and not unbound muons, 
was observed in KC1 (also in NaCl, GaAs). However there is an important difference in the 
nature of the relaxation in diffusing muonium and diffusing muons. In the former case, 
the relaxation rate is given by Tf 1 ~ (^,7^/(1 + col 2 T d) where t^ 1 [27H29] is the diffusion 
rate (inverse 'time-of-stay' of Mu on one site) and S ex is an average electro-nuclear coupling 
constant S ex ~ co c [nI(I + l)/^] 1 / 2 , u c being the nearest-neighbor atoms nuclear hyperfine 
parameters, / the atomic spin, and n the number of nearest neighbors. The energy gap uj\2 
is the smallest intratriplet transition of the Mu spin state related to the (isotropic) contact 
interaction Ui SO by ui 12 = Wj SO [l + (r_/r + )x — (1 + x 2 ) l l 2 }/2. Here T± = (7^ ± 7 e )/2, where 
I11 {le) stands for the muon (electron) gyromagnetic factor, x = 2T + B / Ui SO is a parameter 
which varies with the applied longitudinal field B so that measurements at different fields 
allow one to extract 5 ex and finally the diffusion rate tJ 1 from the raw data [27H30] . The 
value of uj iso were estimated to be about 4280 MHz [33] , while S ex is about 50 MHz in KC1 
(almost temperature independent [28J). Since u iso is large, close to the minimum of t^ 1 
as function of temperature, we have ui2Td S> 1, so that the relaxation rate is found to be 
roughly proportional to the diffusion rate: Tf 1 ~ {^lx/ C0 i2) T d ■ 

By contrast, in the case of a diffusive muon, Kondo [31] found that T{ 1 = co 2 Td indicating 
a motional narrowing decay provided Wo r d ^ 1 (which is the case presently). Thus, for Mu 
the relaxation rate is (roughly) proportional to the diffusion rate, whereas for a muon it 
is inversely proportional, ujq is the second moment of the distribution of fields due to the 
magnetic ions in the compound. In spin ice it is generally considered to be large, and we 
have fixed uq ~ 1 T throughout this paper [TJ [13]. Most importantly, for short times, 
the asymmetry function A(t) within this mechanism is a Gaussian function of time [31] 
A(t) ~ exp(— w 2 t 2 ), so that at the minimum time of measurement, which is about 10~ 8 s, 
with A(t = 10 _8 s) ~ 0.15 already strongly reduced from the fully polarized limit, as is 



observed experimentally. From the relaxation rate data as a function of the temperature 
(Fig.3 in Ref. [I]), we have extracted 1/r^ = 0JqT\ which is the diffusion rate of the muon 
assuming our basic hypothesis. This is represented in our Fig. [2| together with the diffusion 
rate previously measured in KC1 [26T - T28] . We see that the resemblance is quite appealing : 
the diffusion rate is of the same order of magnitude, its range of variation is about 2-3 orders 
of magnitude in both cases, and finally the cross-over temperature (minimum of the curve) 
is more or less the same (50 K for spin ice, 70 K for KC1). An important point to note is 
that the value of tJ 1 at the true minimum (around 50 K) in the case of spin ice cannot be 
measured : this is due to the fact that it would lead to values of Tf which are beyond the 
experimental limit of ~ 10~ 8 s pQ. 

Now, we are left with the problem of understanding the polaronic diffusion in the sub- 
stance as a function of temperature. This problem, which belongs to the general topic of the 
quantum diffusion of heavy particles in solids, has been widely studied in the [311 I37T4J5] . 
All microscopic theories agree that polaronic diffusion occurs in three temperature regimes 
(see Fig. [I]). First, the very low temperature regime (lower than about 7 K in our case) is 
characterized by a diffusion rate independent of temperature: this is the tunneling regime 
of the polaron, where the muon surrounded by its phonon cloud tunnels through the lattice 
as a whole as if it were a rigid free particle. This is a coherent band motion. The high tem- 
perature regime (above T* = 50 K in our case) is a thermally activated one characterized 
by an excitation energy. In this regime the polaron jumps from site to site. Finally in the 
intermediate regime, the motion is still incoherent but its coherence progressively increases 
as the temperature decreases. 

In the spirit of works considering quantum diffusion in KC1 (following Kiefl [27] and 
Kagan and Prokofiev [26]), we have fitted the diffusion rate with the following formula: 

l = ^ r -Eg/T , V ± (IN 

r d ~ VT 1 + (uo/vi) (T/Q D ) a { } 



The first term dominates at high T > T* and corresponds to the usual activated regime of 
the polaron motion [37] 05] . The second term interpolates the very low temperature regime 
with the intermediate one [26J. In the intermediate temperature range, 1/r^ ~ ^{T/Qd) 
[27], whereas at low temperature 1/r^ ~ uq corresponds to the tunneling regime. 

Our fit is shown in Fig. [2j with the values of the different parameters used. The parameters 
appropriate for acoustic phonons give a = 3 (see supplemental material) and Qd = 350 K 
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FIG. 2: Comparison of the diffusion rate (inverse correlation time) t^ of (a) DTO extracted 
from the data of Ref. pQ by using t^ = (VqTi, where we have fixed ujq = 1 Tesla and (b) KCI 
taken from Ref. [26 ! . The original KCI data may be found in Ref. |27| [28] . The fit in panel (a) 
follows the formula Q with a = 3, v = l/r = 1.7 • 10 10 s' 1 , E a = 450 K, v x = 6.8 • 10 11 s -1 K 1 /2 
and vi = 3.5 • 10 5 s _1 . @d = 350 K is the estimated Debye temperature (value taken from |35j). 



as quoted in the literature [32] • This leaves four fitting parameters z/ = 1.7 • 10 10 s _1 , 
v\ = 6.8 ■ 10 11 s~ 1 K 1 / 2 , z/2 = 3.5 ■ 10 5 s _1 and E a = 450 K. By comparison, in case of 
KCI, the different parameters were found to be [27J: E a ~ 390 K, z/ ~ 5 • 10 10 s _1 , u 2 ~ 
1.3-10 6 s _1 . The parameter v\ cannot be directly compared because in [27] a pure exponential 
law Ui exp(—E a /T) was used, while our (ui/\/T) exp(—E a /T) law is more appropriate for a 
polaron. At T = 100 K, we find v x jyfT = 6.8 • lO^s -1 , while the value in KCI is 8.3 ■ 10V 1 . 



The different parameters are roughly of the same order in DTO and KC1. Although Mu is a 
neutral particle, it is is a composite one (muon + electron), whose parts interact differently 
with both acoustic and optical phonons (muons and electrons have quite different masses, 
and thus different effective particle-phonon interaction). We do not discuss this point any 
further, but we think that both types of phonon may also be important for the diffusion of 
Mu in KC1 and that is the origin of the similarity observed between the (possible) diffusion 
of muon in DTO and of Mu in KC1. We note however that to our knowledge, the composite 
nature of the Mu has never been taken into account for its diffusion in any microscopic 
model. 

Equation [I] contains four fitting parameters. These can be related to the parameters of 
a microscopic model (see Supplemental Material for details) (33, HO] to provide estimates 
of their values, which we find to be entirely consistent with the fit we have extracted from 
the experimental data. To summarize the content of the model, we have three distinct 
temperature regimes. At high temperatures, there is a thermally activated hopping regime, 
whereas at very low temperature there is a coherent tunneling regime. Between both, there 
is a third intermediate temperature regime for the hopping rate which exhibits a ~ T -3 law. 

At high temperatures, the hopping rate is: 

' ' ,A 2 e- E ^ kBT . (2) 



Td h^E a (k B T) 

Here E a is the polaron activation energy which is of the order of the typical phonon energy 
[33 SO] • A is the tight-binding hopping term of a muon, related to the total bandwidth 
W ~ ZAq (Z is the coordination number of the lattice). Since typical bandwidth for the 
electron is of the order of a few eV, the bare muonic bandwidth may be estimated to be of 
the order of (m e /m M ) ~ 1/200 times the typical electronic bandwidth. That gives Ao ~ 10 -3 
eV. The coefficient of (1/a/T) exp(— E a /ksT) is just Vy. On the basis of our estimated A 
and the fitted activation energy E a ~ 450 K, we find that v\ ~ 10 11 s -1 K x / 2 , which is of the 
same order of magnitude than the corresponding fitting parameter. 

In the intermediate temperature regime (T < T*), the hopping rate becomes 

I Ag ( t y 3 

r d Hk B Q D ) \e D J ■ { ] 

This fixes the coefficient V2 and a. Then v-i = Aq/^(A;b©c). In this formula, Ao is the 
polaron bandwidth which is strongly reduced from the bare muonic bandwidth by a reduction 



factor exp(— S): A = e~ 5 A . Usually, e~ s is estimated to be 1CT 2 — 1CT 4 (see Ref. [37]). 
By taking 10 -3 , we obtain u 2 ~ 10 5 s _1 in agrement with the fitting parameter. 
Finally, at low temperatures, there is a coherent regime for which 

- ~ A /n. (4) 

Td 

This corresponds to the fitting parameter vq. Since we have already estimated Ao ~ 10 _3 Ao 
we may directly check that it gives uq ~ 10 10 s _1 , also in agreement with its fitted value. 

Our main conclusion is that the measurements of Dunsiger et al. [1] appear to be com- 
pletely compatible with the observation of quantum diffusion of a muonic polaron in its 
three regimes of temperature. Indeed, the magnetism in this particular material apparently 
offers a beautiful probe of the motional crossover phenomenon as the temperature is varied. 
We also propose an experimental check of this basic hypothesis |16] : one may carry out 
the same /iSR experiment on non-magnetic compounds within the same family as DTO 
such as Y2Ti207 and L^T^CV If the mechanism we put forward in this paper is correct, 
muon spin relaxation should also occur due to muon diffusion, but now through the local 
fields of the nuclear spins. (If no nuclear spins were present, muon diffusion would become 
unobservable in /xSR.) If the muons are localized, there would a priori be no possibility of 
relaxation anymore (since the nuclear-spins have very large relaxation times), except by the 
Kubo-Toyabe mechanism which, on its own, would generate an asymptotic 1/3 tail in the 
asymmetry. We have found only one example of such an experiment in the literature (17] , 
in which a relaxation plateau with a relaxation rate of about 1 yus -1 was indeed measured. 
This is consistent with the basic hypothesis that the muon spin relaxation comes from the 
diffusion of the muon. However, we think that further investigations are needed to clarify 
this very interesting physics. 

As a final remark, we note that constant relaxation rates at low temperatures have 
been almost systematically observed in other frustrated magnets. Typically, the constant 
relaxation rate gives way, at higher temperatures, to clear signs of the intrinsic magnetism 
in the material. For example, one often observes signs of magnetic phase transitions using 
/iSR in agreement with other experimental probes. While quantum diffusion may offer an 
explanation for the low temperature relaxation plateaux beyond the case of DTO, there is 
often no sign of the two higher temperature regimes discussed in this paper. However, two 
time scales will be relevant to the experimentally probed relaxation rate: the diffusion rate of 

8 



the polaronic muon Td and 77, the characteristic fluctuation time for the magnetic spins. Both 
quantities vary with temperature. When 7/ <C Td, the muon relaxation is mainly driven by 
the spin fluctuations, and the muon diffusion may be ignored in that case. On the contrary, 
when Td <C Tf, the muon diffusion drives the relaxation. This is probably what happens 
at low temperature in the pyrochlore compounds, which could explain the systematically 
observed plateau at very low temperature. In the spin ice material considered in this article, 
our scenario supposes that the spin dynamics gives a sub dominant contribution to the muon 
relaxation at low temperature and that r^ < Tf for all temperatures. Finally, we note that 
a general theory which includes effects of both r^ and Tf simultaneously seems to be absent 
in the literature [48J, although McMullen and Zaremba [19] partly discussed this case. 

We would like to thank the following physicists for very useful discussions : B. Barbara, 
S. Dunsiger, A. Keren, and T. Uemura. We also thank C. Castelnovo, M. Gingras and S. 
Sondhi for these and collaborations on much related work. 

Appendix A: Supplementary Material 

To represent muon diffusion through the lattice, we start with the following tight-binding 
Hamiltonian [ID] : 

qA n,<5 

1 EE^aV**, (ai) 



N qA n 

The first term is the phonon Hamiltonian and the second term corresponds to the band 
motion of the muon (5 is a vector connecting two neighboring sites of the lattice). Finally, 
the last term is the muon-phonon coupling. Since the muon is a charged particle, it interacts 
with both the optical and acoustic phonons (A is a phonon branch index). For the sake 
of simplicity, we also assume in the following that the muon sites form a simple three- 
dimensional cubic lattice. Let us now introduce the lattice operators: 



M Q,A = K,A + a -q,A)/v2 

P q ,A = (a_ q , A - a+ x )/V2i. (A2) 



As is usual in the theory of polarons, we also introduce the following operators: 

' E— W q - R ", (A3) 



'N ^ uj qX 
= E^nC+c n , (A4) 

n 

and we perform the Lang-Firsov unitary transformation on the Hamiltonian (1) |4UI .'' ], 

e^He-^ (A5) 

The resulting transformed Hamiltonian is: 

H' = ^E^vOpqaP + KaI 2 ) -E a Y, c t c n 

qA n 

- AoE^- M Cc n (A6) 

n,<5 



with 



E a = ^E h f^- (A7) 

2iV qA w qa 

Let us write the last term of the expression (6) in a different form: 

n,<5 n,5 

+ A E [c*<*»-^ - ( e '(*--*-+*)> J c + +5 c n . (A8) 

The expression (^4) , means that we take the mean value of the operator A on the phonon 
states | . . . A r kA . . . N k >\ . . .), where N k \ are the phonon occupation numbers, 

^ = ps^n < A9 > 

The expression (8) defines two different channels for the polaron motion from one site to 
the neighboring site: a coherent channel where the phonons occupations numbers do not 
change during the motion (diagonal term in (8)), and an incoherent channel (jump motion) 
where changes in the phonon occupations numbers during the motion are allowed. At very 
low temperatures, the second process can be neglected, whereas it becomes predominant at 
higher temperatures. The diagonal term can be calculated and gives: 

(e*(*»-*»+<)) = e s(T) ( A10 ) 
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S(T) = ^ E ^# ™th (/3u;qA /2) sin 2 [(q • 5) /2]. 

7V qA W qA 



(All) 



Coherent motion. Let us first examine the coherent motion at very low temperatures. In 
this regime, all non-diagonal phonon transitions are negligible so that the relevant effective 
Hamiltonian for the polaron motion can be reduced to: 



H, 



coh 



-£a£c+c n -A oe -^>£c+ +5 c n . 



(A12) 



n,<5 



The eigenstates are e(k) = — E a — 2A [cos(k x ) + cos(k y ) H-cos(fe^)], with A = A e S ( T > 
and the Green's function is 

(A13) 



G q (;) = ^E ei(ek ~ ek+q)< - 



From the above expression we can calculate the mean time of stay of the muon on one site, 
i.e. its inverse tunneling rate: 



r 






r (E ei(£k ~ £k+q)4 ) dt 



h h 

~ — rsJ — 

2V/2AQ A 



(A14) 



o 



The expression (11) for S thus defines the polaron reduction factor of the bandwidth in the 
coherent regime. 



Incoherent hopping rate: At intermediate and large temperatures, the phonon occupation 
numbers Nf. increase so that the second channel for the diffusion becomes predominant. At 
these temperatures, the delocalized states k lose their meaning and the effective Hamiltonian 
for the polaron diffusion is reduced to: 



H incoh = -E a E c+c n - A E [ e * ( 



K+S) _ I p l(4>rx-<t>n+6) 



n.d 



ph 



C n+<5 C n- 



(A15) 



From the second term, we can calculate a transition probability W(n — > n ± 5) to second 
order in perturbation theory, and the corresponding jump (or hopping) rate. Holstein carried 
out this calculation and got the following hopping rate [37] : 



1 
r 



A 2 

1/2 

7T ' 



h 

x exp 



E |7ka| 2 sin 2 [(k ■ 8)/2] csch (/3^ kA /2) 



kA 



-1/2 



- 2 E T^" sin2 [(k ■ S)/2) tanh (/3w kA /4) 



kA W kA 



(A16) 



11 



(we recall that csch(x) = sinh(x) -1 ). 

In the high temperature limit, i.e. T ^> u/2, we have tanh (/3o>/4) ~ 0uj/4 and 
csch {j3uj/2) ~ 2//3oj, and using the definition (7) of E a , (16) becomes 

1 7T 1 / 2 

_ A 2 e -E a /k B T^ (A1?) 



t h^/E a (k B T) 

which corresponds to our expression (4) in the paper. This is the thermally activated regime. 
In the intermediate temperature regime which should nevertheless be a low temperature 
with respect to the phonon frequencies, tanh(/3co>) ~ 1 and csch(/3co>) ~ exp(— f3uj) so that 
(16) becomes by using the definition (11), 



1 A 2 

r h 



$> k A| 2 sm 2 [(k-5)/2]e- 

kA 



-(/3a, kA /2) 



-1/2 

(A18) 



To evaluate this expression, we must now separate the respective role of the optical and 
acoustical branches. The optical branches give the main contribution to the reduction factor 
S(T) [37], whereas the contribution to this factor coming from the acoustical branches is 
much less J43] . However in this intermediate temperature regime, the acoustic branches play 
a major role in the polaron motion. The reason is that the exponential factors in (18) vanish 
for all optical phonon frequencies. Consequently, we may only consider the acoustic phonons 
in the summation over the different branches in (18). Let us now calculate the factor 

£ | 7kac , I 2 sin 2 [(k acc . • 8) /2\ e -0**~/*> (A19) 

K-acc. 

We adopt the Debye approximation for the acoustic phonons, and following Petzinger [40J, 
we take for all longitudinal modes 

Eacc. w |7k| 2 /2w k , 
sin 2 [(k • 8)/ 2} w luI/lu 2 d , (A20) 

where ujd is the Debye frequency (Tiujd = ^b@d = koc/h, with 6^ the Debye temperature, 
c the sound velocity and ko ~ 7r /a the Debye wave vector). By introducing the phonon 
density of states in the Debye approximation go(w) = [3/(27r 2 )] u 2 /uf), and replacing the 
summation in (19) by an integral, we easily get 

*~ , A ° Gfr. (A2D 

12 



where all the numerical factors together have been evaluated to be of the order of unity. 
The precise value of -E acc . is not known, but is necessarily of the order of fcg0£), so that we 
have finally 

l~W^)(£) ' (A22) 

which is the formula (3) of our paper, which is applied in the intermediate regime. 

Finally, we have adopted the simplest theory of the quantum diffusion of a muonic polaron 
to justify our fit, and it appears quite satisfying. However, other possibilities such as two- 
phonon interaction processes and/or a fluctuation preparation barrier [3T] should be also 
examined. Possible Umklapp processes, when several muon sites per unit cell exist, should 
also be studied. We reserve all these technical discussions for a further publication. 
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